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Abstract 

The subject of the first section-lecture is concerned with the strength and the 
weakness of the perturbation theory (PT) approach, that is expansion in powers 
of a small parameter a, in Quantum Theory. We start with outlining a general 
troublesome feature of the main quantum theory instrument, the perturbation ex- 
pansion method. The striking issue is that perturbation series in powers of a <^ 1 
is not a convergent series. The formal reason is an essential singularity of quantum 
amplitude (matrix element) C(a) at the origin a = . In many physically impor- 
tant cases one needs some alternative means of theoretical analysis. In particular, 
this refers to perturbative Quantum Chromodynamics (pQCD) in the low-energy 
domain. 

In the second section-lecture, we discuss the approach of Analytic Perturba- 
tion Theory (APT). We start with a short historic preamble and then discuss how 
combining the Dispersion Relation with the Renormalization Group (RG) tech- 
niques yields the APT with e~^^" nonanalyticity. Next we consider the results of 
APT applications to low-energy QCD processes and show that in this approach 
the fourth-loop contributions, which appear to be on the asymptotic border in 
the pQCD approach, are of the order of a few per mil. Then we note that us- 
ing the RG in QCD dictates the need to use the Fractional APT (FAPT) and 
describe its basic ingredients. As an example of the FAPT application in QCD 
we consider the pion form factor F.^{Q^) calculation. At the end, we discuss the 
resummation of nonpower series in (F)APT with application to the estimation of 
the Higgs-boson-decay width r^^^^(m|^). 
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1. Strength and weakness of perturbative QFT 

1.1. Essential singularity at a = and functional integrals 

In quantum mechanics and QFT we have a lot of successful perturbative 
calculations. Practically, Perturbation Theory is a synonym of Quantum 
Theory. Feynman diagrams became a symbol of QFT. Nevertheless, pertur- 
bative power expansion of the quantum amplitude C{a) is not convergent. 



Feynman Series ^ c^a'^ is not Convergent ! 

The reasons for this behavior are rather simple. First, remind the Dyson 
argument pT] . In QED the substitution a — >• —a is equivalenlQ to e — t- z e 
and due to 

this change destroys hermiticity of Lagrangian and unitarity of S-matrix. 
Hence, in the complex a plane, the origin cannot be a regular point, instead, 
one has an essential singularity at a = 0. 

Second, in QFT one meets factorial growth of coefficients ~ kl and 
this is due to an ill-posed problem. A small parameter g standing at the 
highest nonlinearity is an indispensable attribute of quantum perturbation. 
Indeed, we quantize a linear system as a set of oscillators. Only after 
that we account nonlinear terms ~ 5 <C 1 as small perturbations. However, 
nonlinearity usually changes the equation seriously — new solutions appear. 

The most general way to analyze the issue is to use the functional 
integral representation. Here, for illustration purposes only, we consider 
the analog (often called "0-dimensional" ) of the scalar field theory gip'^ 



I{9) 



00 



-gx^ 



dx . 



(2a) 



Expanding it in power series 

T, ^ ^kr -.u T r(2A; + l/2) 

Ii9) = Y.{-g)'h with 4= r(fc + l) 



k=0 



2^ k\ , (2b) 



fe>i 



one arrives [2] at factorially growing coefficients. Meanwhile, I{g) can be 
expressed via special MacDonald function 



(2c) 



with known analytic properties in the complex g plane: It is a four-sheeted 
function analytical in the whole complex plane besides the cut from the 



^Here a = e^/(47r) is the QED expansion parameter, the fine structure constant. 
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origin g = along the whole negative semiaxes. At the origin, it has an 

essential singularity e~^^^^ and in its vicinity on the first Riemann sheet it 
can be written down in the Cauchy integral form: 

(2d) 




7(5 + 7) 

As far as the origin is not an analytical point, the power Taylor series ()2bp 
has no convergence domain for real positive g values. This is in concert with 
factorial growth of power expansion coefficients. The power series is not 
valid also for negative g values — in accordance with Dyson's reasoning. 

Besides, via 0-dimensional analog of functional integral (j2ap one can il- 
lustrate the analysis of the /(g() analytic properties in the complex g plane 
by the steepest-descent method which was devisecQ for the functional in- 
tegral representation. Then it is possible to prove [3] factorial growth of 

expansion coefficients in the (/>^ scalar and a few other QFT models. These 
results have been anticipated in 1952-53 [4| just after Dyson's paper. 

The same singularity structure ~ exp(— was established by im- 
proving the perturbative result using two other nonperturbative methods: 
Q^-analyticity (from Dispersion Relations), and Renormalization Invari- 
ance [5j: 

/pert(Q',5) = l-/3o5ln(Q') ^ /imp(Q'e-^/'^o^'). (3) 

Henry Poincare analysis of Asymptotic Series (AS) properties at the 
end of the XIX century can be summarized as follows: 

/j. The truncated AS can be used 

for obtaining quantitative informa- 
tion on expanded function. Here, 
^ the error of approximating F{^g) by 

first K terms of expansion, F{(j) — > 
• FK{g) = ^k<K fkig) is equal to the 

last detained term fxig)- 
• 9 This yields to the existence of a 

( ^ lower limit of possible accuracy for 

the given g value (in contrast to 
1 2 3 ... K K + i k convergent series!). 



To elucidate the phenomenon, take a power AS, fk{g) = fkg^ with 
factorial growth ~ /c! of expansion terms fkig) which cease to diminish 
at k = K ^ 1/5 - For k > K + 1 truncation error starts to grow! 

This yields the natural best possible accuracy of an AS at a given value 
of expansion parameter. From the explicit illustration for the function I[g) 
(|2ap with AS ()2bp . presented in Table [H one can see that the optimum 
values of truncation number K = ~ l/{2g) (in dark blue) provides us 



^See, e.g., Section 2 in the paper [2] and references therein. 



4 



Alexander P. Bakulev, Dmitry V. Shirkov 



with the best possible accuracy (dark blue in the last column). Indeed, an 
attempt to account a couple of extra terms (the second and the forth lines) 
results in drastic rise of the error! 



.9 


K 


(-.9)^/k 


(-g)^+i/K+i 


lK{g) 




Ak/(5) 


0.07 


7 


-0.04(2%) 


+0.07(4.4%) 


1.674 


1.698 


1.4% 


0.07 


9 


-0.17(10%) 


+0.42(25%) 


1.582 


1.698 


7% 


0.15 


2 


+0.13(8%) 


-0.16(10%) 


1.704 


1.639 


4% 


0.15 


4 


+0.30(18%) 


-0.72(44%) 


1.838 


1.639 


12% 



Table 1: The last detained {(— g)^ Ik) and the first dismissed {{—g)^^^lK+i) 
contributions to the power series (j2b|) in comparison with the exact value (|2cp and 
the approximate result Iidg), with K being the truncation number and Afclig) 
— the error of approximation. 

Thus, one has i^*((? = 0.07) = 7 and K^{g = 0.15) = 2 . It is not possible 
at all to get the 1% accuracy for g = 0.15 . 

In QED this 'divergence menace' is not actual, as the real expansion 
parameter is quite small: a/vr ~ l/(1377r) ~ 2 • 10""^ . At the same time, in 
perturbative QCD (pQCD) the expansion parameter below 5-10 GeV is not 

very small: as{Q) ~ 0.2 — 0.3. Then for an observable Aqcd = J2k O'k{os)'' 
with ttk ~ k\ and with critical order ~ 3 — 5 the danger of the pQCD 
series explosion is actual, see below in Tabled 

Hence, a practical Non-Perturbative approach is of utmost importance. 
In Section 2, we concentrate on the so-called Analytic Perturbation Theory 
(APT), a closed scheme, devised in the late 90s[B]. It combines informa- 
tion from PT with two other nonperturbative methods — Analyticity and 
Renormalization Group (RG). Having this in mind, we outline now the RG 
approach. Then ideas of the Dispersion Relation method will be presented. 

1.2. RG transformation 

Consider transformation Rt [in — t- fik ■, di ~^ Ok] as operation with continu- 
ous positive parameter t, acting on a group element Qi = Q{fii,gi), specified 
by 2 coordinates fii and g^. This operation 

RfQi = Qk^ Rt {^i /ifc = tfii, gi^ gk= g{t, gi)} (4) 

contains dilatation of fi and functional transformation of g'^ . The Rt group 
structure Rt Rt = Rtr is provided by the functional equation (FEq) 

9{rt,g) =g{T,g{t,g)) . (5) 

Indeed, if one puts x = rt, then its LHS describes Rrt acting on g, Rrt 9 = 

'g{Tt,g) , while the RHS one corresponds to the two-step procedure: 

Rr ® Rtg = Rt g{t, g) =g (t, g{t, g)) . We see that a combination of the two 



'^For a more detailed exposition of tliis material kindly address to [?]■ 
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lines results in Eq. (jS]) , providing the group composition law Rrt = Rt "X" Rt ■ 
Thus, the operation Rt forms a continuous Sophus Lie (1880) group of 
transformations. 

The RG symmetry and RG transformation are close to the notion of 
self-similarity, well-known in Mathematical Physics since the end of the 
XIX century. The Self-Similarity Transformation (SST) is a simultaneous 
power scaling of arguments z = {x, t, . . .} and functions Vi{x, t, . . .) 

Sx:{x^ Xx,t^ \U} ; {Vi{z) ^ Vl{z') = X'^Viiz')} (6) 

Below, we call it the Power Self- Similarity (PSS) transformation. 

The general solution of g{xt,g) = 'g(x,g{t, g)) depends on an arbitrary 
one-argument function see below Eq. dS]). Here, we look for a partial solu- 
tion, linear in the second argument 'g{x,g) = g ■ f{x). The function f{x) 
satisfies simple FEq f{xt) = f{x) ■ f{t) with general solution: f{x) = x^ 
and ^(x, t) = g ■ x^ . Thus, RG transformation is reduced to the PSS one. 



Rt ^ {x ^ x - t, g ^ g ■ t"} = St. 



The PSS transformation Rt ^ St = {x ^ x ■ t, g ^ g ■ f^} is a special 
case of the RG one. That is, in the RG case instead of the power law , 
one has arbitrary functional dependence. Hence, one can consider the RG 
transformation as a functional generalization of the PSS one. It is natural 
then to treat them as transformations of functional scaling or Functional 
Self- Similarity (ESS) transformation. In short RG = ESS. 
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Here we illustrate our statements by historical analogies. One realizes 
that scaling = change of scale = proportional change of sizes. This notion 
implies that Bogoliubov Renormalization Group (BRG) = distorted scaling 
= change of scale with continuous change of some details. This is illustrated 
by Eig.[Tl in both paintings one can recognize the same person, Leonard 
Euler but details are differen10. 



''The original (smaller in Fig.[T| portrait was painted by E. Handmann in 1756 and 
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The symmetry of the FSS group transformations can be 'discovered' 
in different fields of physics. As the illustration, we suggest a mechanical 
example. Imagine an elastic rod with a fixed point (point "0") bent by 
some external force, e.g., gravity or pressure of a moving gas or liquid, see 
in the left panel of Fig.[2j The form of the rod can be described by the 




Figure 2: Elastic rod model. 



angle g between tangent to the rod and vertical directions considered as a 
function of the distance / along the rod from the fixation point — by the 
function g{l)- If the properties of the rod material and exterior forces are 
homogeneous along its length (independent of/), then g{l) can be expressed 
as function G{l,go), depending also on go — deviation angle at the fixation 
point from which distance I is measured. G can depend on other arguments, 
like extra forces and rod material parameters, but in this context they are 
irrelevant. 

Take two arbitrary points on the rod, "1" and "2" with h = X and 
l2 = X + l. The angles gi at points "0", "1" and "2" are related via G 
function: 

gi = G(A, go), 52 = G{X + I, go) = G{1, gi) . 

To get the RHS of the second equation, one has to imagine that the fixation 
point now is "1", as shown in the right panel of Fig.[2j Combining both 
the equations, gi = G{X,go) and g2 = G{X + I, go) = G{l,gi), one gets the 
group composition law G{X,G{l,g)} = G{X + l,g) equivalent [via relation 
'g{x,g) = G{1 = Inx, g)] to FEq for invariant coupling 

9{x,g)=g{j,g{t,g)). (7a) 

first devised in 0. It can be represented in two infinitesimal forms: The 
nonlinear differential equation (DEq) 

m=t^ (7b) 



dx dt 



t=i 



is stored at Basel University, whereas the larger and more official one was painted by 
I. Konig at the request of the Russian Academy of Sciences (due to its 150-year anniver- 
sary) in 1875 as a copy from the original portrait. This newer portrait was lost during 
Russian revolution and discovered in a curiosity shop in 1972 by Georgy Sergeevich 
Golitsyn [8]. 
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and the linear partial one (PDEq) 



dx ^^^^ do 



g{x,g) =0. 



(7c) 



In the QFT jargon, the beta-function [5{g) is known as the RG generator. 
Note that each of the two DEqs is equivalent to the functional one (I7a|) 
provided the normalization condition 



5(1,5) = a 



(7d) 



is satisfied. The general solution of the last FEq can be found from the 
relation 



J a 



3 dj 



Inx . 



(8) 



1.3. RG efTective coupling in QFT 

To illustrate the power of RG-invariance, on the one hand, and the effective- 
ness of its differential formulation, consider now an illuminating example 
of the effective coupling g UV asymptotics. In the common perturbation 
theory (PT), one has 

9pt{x; g) = g + g'^f^o In x , (9a) 



the so called one-loop UV logarithm. Here, as well as in Eqs. ([7al 
x = Qyfi^- g = g^ = g{x-g); Q^ = Q^-Ql>0. 



Evidently, expression (|9a|) is not RG-invariant. Indeed, substituting it 
into the functional equation ([5]) one obtains the discrepancy 



A 



discr 



5pT 



5pt(^; 9) - 9pT (f , 9pt(*; 9) 

[g + g^Po Inx] - [g + g^(3o Inx + 2g^(3l Int \n{x/t)] / 



— error of the (7^-order that can be killed by adding to starting approx- 
imation (|9a|) the next-order term g^ /3q In^ x etc. The final result of this 
iterative restoring is the famous sum of geometric progression 



fc>0 



Inx 



(9b) 



which sums up the leading order (LO) logarithms g{g Inx)*^ . On the other 
hand, this solution can be immediately obtained by analytical means via 
the first differential RG Eq. ()7bp with /3 (g) = /3o 5^ > obtained from Eq. ([9a|) . 
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Starting with the two-loop perturbative UV asymptotics 

g)=g + g^Po In x + (^32 (^^^ ^ ^-^q^^ 

one gets (again in few hnes of calculations by RG technique) the Next-to- 
Leading-Order (NLO) approximate (at Inx ^ 1) result 

^Rg(^;5) ^ 1 ^1 ^ 2, , n , ■ (10b) 

1 — g Po Inx — g"^ Pi ln(ln x) 

Illustration in QED. In QED, due to gauge invariance (Ward identi- 
ties), the RG-invariant coupling reduces (see, pioneer review paper [9] and 
references therein as well as Sect. 48.1 in monograph [TU]) to the transverse 
amplitude of the dressed photon propagator 

a{Q'^,a) = adtr{Q'^,a) . (11) 

Originally, it was introduced there as a function describing the RG trans- 
formation of a PT expansion parameter. In a proper QED (quantum elec- 
trodynamics of electrons, positrons and photons) its one-loop expression 
reads ^2 

apUx;aa) = au. + ^ Inx; x = ln((5^/^^) . (12) 




°^0> 7.137 
. Millikan 



100 GeV 



Figure 3: Artistic view on experimental verification of the oqed "running" (solid 
line) . 



In a sense, it is a generalization of the electron effective charg^EI intro- 
duced first by Dirac [TT] : 

„2 



e{Q) ~ e 



1 + 



^ In 

12 vr^ Xnie 



+ ... 



After the RG machinery, ()12p takes the invariant form 



aQED(Q^ 



a. 



l-(a^/(37r)) ln(QV/x2)- 



(13a) 



(13b) 



^Its Fourier image is the charge Q{r) of point electron screened by quantum vacuum 
fluctuations. 
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Its contemporary Standard Model analog was checked experimentally about 
a decade ago at LEP, the e"'"e~ -collider at CERN, see results in Fig.[3l 

Note also that the Bogoliubov RG was initially devised in a mass- 
dependent form [3. Starting, e.g., with the mass-dependent one- loop PT 

g^U^, y;g)=9 + g^ \h {x/y) - h (i/y)] (I4a) 

with X = / fi^ and y = m^//i^ , one gets 

^- g [h{x/y) - Ii{l/y)\ 

and, similarly at the two-loop level|12j. This massive RG provides us with 
means to make an accurate matching across the heavy-quark thresholds, 
both in QCD and in Grand Unification. 

1.4. Analyticity from causality in QFT 

Turn now to the Dispersion Relation Method that relates causality in space- 
time with analyticity in kinematic (energy, momentum transfer) variables. 
To illustrate the main idea, consider the Fourier image 

00 

F{E)= [ e''^A{t)dt (15) 



of the forward scattering amplitude A{t), being subdued to the nonrelati- 
vistic causality condition: 

A{t) = at t<0. (16) 

In this case, F(E) can be analytically continued from real E values to the 
upper half of the complex plane Q ^ z = E + i^] ^ = Imz > 0, since the 
factor e~*^ in the integrand provides convergence of the integral for F{z). 
Using then the Cauchy theorem with integration contour F"*" in the upper 
half-plane and z inside F^ 

ZTT I Jy+ Z' — Z 

one can get Dispersion Relation for the forward scattering amplitude 

00 00 

—00 m 

relating two observable functions. In obtaining this non-subtracted dis- 
persion relation we tacitly assumed "good" F{z) asymptotic behavior and 
used Optical Theorem ImF{E) = ka{E). In a more realistic case, one 
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starts with relativistic causality and adds symmetry crossing property of 
the forward scattering amphtude. 

Remind that instead of the Pauh-Jordan commutator 

Z)(x-y) = i (0|[(/.(x), (/>(y)]|0), (18a) 

(vanishing outside the hght cone at (x — y)^ = {xq — yof' — (x-y)^ < 
and involved in relativistic-invariant quantization) in construction of ma- 
trix elements and observables one needs the Stueckelberg-Feynman causal 
propagator 

D,{x -y) = D^{x-y) = -^ \T [(/.(x) 0(y)] | 0) , (18b) 

being the vacuum expectation value of the time-ordered product. Just this 
function (and its derivatives) enters into Feynman rules. 
At the same time, in the Schwinger-Dyson equations one deals with its 
"dressed" version which includes radiative corrections 

D^rU^; 5) = -V (0 1^ ['^(^) '^(0) S{g)] I 0) , (18c) 

from the scattering matrix 

S{g) = l + gSi + g^S2 + ... 

This dressed causal propagator can be represented in the Kallen-Lehmann 
spectral form0 

oo 

Ddvcssiq ,g) = 2 — ~ + ~ / 2 — ~ (^^"i) 

— q'^ — le vr J a — q"^ — \£ 

mi 

This representation, in a sense, resembles the forward dispersion relation 
(fT7|l . Here, quite similarly, it defines the function D{z) = D^rcssiz, g) 
analytic in the whole complex plane z except the pole and the cut on part 
of the real axis. In its proof, as well, the relativistic generalization of the 
causality condition (fT6]l is used. 

In what follows, we have to combine two nonperturbative methods — 
the RG and the Dispersion Relation ones. 

1.5. RG and causality 

Return to the QED invariant coupling. According to Eq. (jlip . it is propor- 
tional to the transverse amplitude of the dressed photon propagator. As it 

®Here we give its simplest, nonsubtracted version that is valid for the appropriately 
decreasing spectral density behaving as p{a,g) < 1/ln^ cr). 
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was first shown by Kallen f l31, t he latter amplitude can be represented in 
the form congeneric to Eq. (flSdj) . Hence, the same presentation 

oo 

a«3^«)=i /47#^ (19) 

TT y a + -le 



has to be valid for the QED effective couphng. That is, the function a(z, a) 
should be analytic in the duly cut complex z plane. This is true, term- 
by-term, for its PT expansion, li ke E g. (112p . However, it is not true for 
its RG-invariant counterpart, Eq. ()13bp . as it contains parasitic singularity 
outside the allowed cut — the so-called ghost or Landau pole at 

Q2 = Q2^^2g3V«M. (20) 

An elegant solution to resolve the issue was proposed by Bogoliubov et 
al. jl4] . Omitting the details we mention here the main recipe: 

To bring a RG-invariant but singular (containing extra pole) 
expression in the proper Q'^ -analytic form, one has to use the 
spectral representation U9\) wi th th e spectral density p defined 
from the related PT input, Eq. ^3b\) . 

The resulting analyticized expression is 

The second term in the RHS is invisible in the PT expansion. Remarkably, 
it contains an essential singularity at a = of a proper type. 

1.6. QCD effective coupling 

In QCD, at the one-loop, LO approximation one has 

1 + as(^)/3o ln(QV//2) /3qLq J 

The NLO, or two- loop expression in the De nom inator Representation con- 
tains log-of-log dependence (just like in eq. ljlObp for Lq » 1): 

. , _^(o\o,, r ; /3o,i~l. (23) 
^qLq + (/3i//3o)lnLQ 

Note that the famous Asymptotic Freedom UV behavior as{Q) ~ is 
correct already in the LO approximation. The QCD scale A turns out to 
be numerically close to the confinement scale A ~ 300 — 400 MeV ~ 2 771,^; 
that is Ra ~ 10~^^ cm. 
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The QCD final product is formulae for hadronic observables. Some of 
them, e.g., the ratio of inclusive cross-sections for e~^e~ annihilation, can 
be directly expressed 

n / \ '-'^e+e" — ^hadrons \S/ r, / \ /r,A \ 

Re+e-{s) = -— = R{s]as) (24a) 

in terms of the QCD notions. As an RG-invariant this ratio should depend 
on the QCD coupling as{s) onlyl Perturbatively, it is power functional 
expansion 

R{s; Us) = Rinv{as{s)) = l + n as{s) + ra ai{s) + O (a^(s)) . (24b) 

Remarkably that, according to the 2007 Bethke review [15], above a few 
GeV the two-loop pQCD nicely correlates all as{Q'^) data, see Fig.Hl 




Figure 4: Triumph of the NLO pQCD from 200 up to 3-5 GcV (taken from p]V 

However, data below 5 GeV are not good enough: One can see from 
Table [2] that while the situation with higher-loop corrections for the R^+e- 
ratio is rather good (the first line), the Bjorken Sum Rule (SR) case (the 
third line) delivers to us a signal of possible blow-up (like in Table [l]) of the 
related asymptotic perturbative series at Q < 2 — 3 GeV. Just in this region 



Observable 


Scale 


1-loop 


2-loop 


3-loop 


4-loop 


^cxp 


-^e+e~ — ^-hadrons 


10 GeV 


92% 


7.6% 


1.0% 


-0.6% 


12-30% 


Rr in r-decay 


2 GeV 


51% 


27% 


14% 


8% 


5% 


Bjorken SR 


2 GeV 


56% 


21% 


12% 


11% 


6% 



Table 2: Relative size of 1-, 2-, 3-, and 4-loop contributions to observables. Two 
last lines are taken from [TB]. Three- loop estimations can be found in |17| . 
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the PDG canonized explicit expression for as{Q'^) starts to grow sharply 
as even the one-loop expression for effective QCD coupling 

has unphysical /Landau) singularity. Just this singularity at \Q\ = A ~ 350 
MeV prevents from analyzing data by pQCD in the low-energy physical 

region below few GeV. 

Meanwhile, all nonperturbative lattice-QCD simulations testify to reg- 
ularity of as{Q) behavior in the ~ region — see fresh overviews I18j . 
In what follows, we shall deal with ghost-free analytic QCD couplings (and 
their powers) taken from the APT and its generalization which are regu- 
lar in the low-energy region and in the high-energy one coincide with the 
common as(Q^). 



2. Analytic Perturbative Theory (APT) in QCD 
2.1. APT: Historic Preamble 

As it has been mentioned above in Sect. I1.5.|, the analytization recipe is a 
common product of two nonperturbative methods — the RG and DR ones. 
Originally, it was formulated [13] for QED in the Euclidean region. Then 
in 1982 Radyushkin and Krasnikov&Pivovarov [TO] using the dispersion 

technique suggested regular (for s > A^) QCD running couplings in the 
Minkowski region, namely it~^ arctan(7r/L<j) and <I>(q"[L<,]), = ln(s/A^), 
with indication of the need to use the nonpower expansion in <I>(a"[Ls]) 
instead of the power one. The real proliferation of this technique into QCD 
was initiated by Igor Solovtsov and his co-authors in the mid-90s [6]. In 
these pioneer papers, the ghost-free expressions for the RG-invariant QCD 
couplings in both the energy-like Minkowskian region and the momentum- 
transfer Euclidean one were obtained. Quite soon this construction was 
developed [20] in a closed schemeQ See the history details in the recent 
review paper |17] . The whole construction is known since then as the 
Analytic Perturbation Theory|3 

The next step, made by Bakulev, Mikhailov, and Stefanis, generalizes 
the APT by including fractional powers of coupling, as well as products of 
coupling powers and logarithms ^3] and for this reason, it was named the 
Fractional APT. At the same time, it appears possible to sum up nonpower 
series in the (F)APT [M1E5]. 



2.2. Basics of APT 

In the standard pQCD, as we described in Sect. 1, the one- loop RG equation 
for the effective coupling as{Q^) = a[L\/ f3f with L = ln(Q^/A^) and /?/ = 

^The essential point was discovering the necessity of nonpower type of functional 
expansions [21) which are the only compatible ones with linear integral transformations 
relating Euclidean, Minkowski and space position pictures. 

* Close results were partly obtained by Simonov using the background perturbation 
theory [22]. 
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Figure 5: Contours of integration on the road to the Minkowski space. 



(3o{Nf) = (11 - 2iV//3)/(47r) generates the pole singularity, a[L] = l/L. 

Due to this pole, in pQCD the problem arises: How to go to the 
Minkowski region? Quantities in the Minkowski region are usually rep- 
resented by contour integrals of the type ^ f{z)D{z)dz. In the integrands 
one uses D{z) = X^j„dma™(-z) and changes the integration contour to Fj. 
This change of the integration contour is legitimate if D{z)f{z) is ana- 
lytic inside the circle everywhere. But as{z) and hence D{z)f{z) have the 
Landau pole singularity just inside! In the APT effective couplings A n(z) 
are analytic functions and this problem does not appear at all! In [2S] 
(2010) the equivalence of the Contour-Improved-PT (GIFT) approach to 
the (F)AFT for quantities like R{s) was proved, which can be symbolically 
expressed as 



CIPT 



D{z)dz 



Fa 



(F)APT 



D{z)dz 



Fa 



By the analytization in the APT for an observable f{Q'^) we mean the 
Kallen-Lehmann representation (jl8d|) 



■ da 



(26) 



with Pf{a) = — Im[/(— o")]. Then in the one-loop approximation Pi{a) 
vr 



+ anc0 

ML] 

MLs] 



a + Q^ 



da 



da 



1 

I 



— arccos 
vr 



(27a) 
(27b) 



We use the notation f{Q) and f[L] in order to specify the arguments we mean — 
squared momentum or its logarithm L = ln(Q^/A^), that is f[L] = /(A^ • e^) and 
A^ is usually referred to the Nf = 3 region. Note that here we introduced the notation 
.4„[L] and 2tn[L] in order to distinguish analytic images of normalized coupling powers, 
a" (Q^), from the corresponding images of a"(Q^) — for the latters we use the standard 
notation ^n[i] and 2ln[i]. This notation is different from that used in [231 125] . 
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whereas analytic images of the higher powers (n > 2, n G N) are: 



ML] 



in 



1)1 



d ■ 


n— 1 ^ 







(27c) 



Note that at L ^ 1 the pole remover ~ ~ e~^^'^. In other words, 
Kallen-Lehmann analyticity in the plane generates nonperturbative 
g-i/os correction. This correction guarantees the absence of spurious Lan- 
dau-pole singularity and ensures the correspondence with PT as{Q'^) at 
Q2 > 1 GeV^. 
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Figure 6: Distorting Mirror for analytic couplings in the Minkowski and Euclidean 
regions. Left panel: for 2ti(s) and Ai{Q^). Right panel: for 2l2(s) and A2{Q'^)- 



In Fig.[6l we show the so-called Distorting Mirror for analytic couplings 
in the Minkowski and Euclidean regions: in the left panel — for 2li(s) 
and Ai{Q'^), whereas in the right panel — for 2t2(s) and A2{Q'^)- We see 
that in the IR domain one has universal finite IR values ^i(O) = 2ti(0) = 1. 
Moreover, starting from the two-loop level analytic couplings reveal loop 
stabilization of IR behavior. This yields practical loop- and renormaliza- 
tion-scheme-independence of Ai{Q'^), 2li(s), and higher expansion func- 
tions, for details see [T7] . 

In Fig. [71 we sho w th e results of the APT application in meson spec- 
troscopy obtained in [26j . We observe from this comparison that the three- 
loop Us, shown as dot-dashed line, is completely excluded, whereas the 
three-loop APT coupling goes through all the "experimental" points due 
to well-established P- and S'-wave meson states. 

In the APT, in addition to the regular behavior of couplings, one has 
to have nonpower expansions for physical observables. Indeed, if in the 
standard pQCD for an observable D we hav43 

^pt(Q^) = dQ + di as{Q^) + d2 aliQ^) + d^ al{Q^) + . . . ; (28a) 
Rpt{s) = do + di as{s) + d2 a^(s) + a^(s) + ... , (28b) 

^"Here ra and higher-order coefficients r„ differ from ds and dn by the vr^ terms. 
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Figure 7: Comparing af^^ from tlie meson spectrum and tlie tiiree-loop Ai at 
A[f^^^3 = (417±42)MeV (solid lines). The three-loop as is shown as the dot- 
dashed line. 

then in the APT we should use the nonpower functional expansion 

-DAPTiQ^) = do + di AiiQ^) + d2 A2{Q^) + d3 A^{Q^) + . . . ; (28c) 
^APT (s) = do + c^i 2ti (s) + 2I2 (s) + da 2t3 (s) + . . . . (28d) 

This provides 

• Better loop convergence (compare Tables [2] and [3]) and practical renor- 
malization-scheme independence of observables; 

• Third terms in ()28cp and ()28dp contribute less than 5%, cf. Table [3l 
Again the two-loop (N^LO) level is sufficient. 



Observable 


Scale 


1-loop 


2-loop 


3-loop 


4-loop 


^cxp 


-^e+e~ — ^hadrons 


10 GeV 


92.2% 


7% 


0.7% 


0.1% 


12-30% 


Rr in r-decay 


2 GeV 


90.6% 


8.2% 


1% 


0.2% 


5% 


Bjorken SR 


2 GeV 


75% 


20.5% 


4.55% 


-0.05% 


6% 



Table 3: Relative size of 1-, 2-, 3-, and 4-loop contributions to observables in the 
APT (for comparison with PT — see Tablc[l]in Section [1.6 -I) . Estimates on Rt are 
taken from |27] . whereas those on Bjorken SR — from the paper in preparation 
[l6] (the 3-loop estimations can be found in [T7]). 



2.3. From APT to FAPT 

At first glance, the APT is a complete theory providing tools to produce 
an analytic answer for any perturbative series in QCD. However, in 2001 
Karanikas and Stefanis [28] suggested the principle of analytization "as 
a whole" in the plane for hadronic observables, calculated perturba- 
tively. More precisely, they proposed the analytization recipe for terms like 

JqcIx f^dyas [Q'^xy) f{x)f{y) which can be treated as an effective account 
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Figure 8: Comparing 21^ [i] (left panel) and A,y[L] (right panel) vs. L for frac- 
tional G [2, 3]. 



for the logarithmic terms in the next-to-leading-order approximation of the 
pQCD. Indeed, in the standard pQCD one also has: 

(i) the factorization procedure in QCD that gives rise to the appearance of 
logarithmic factors of the type: a'^[L] L; 

(ii) the RG evolution that generates evolution factors of the type: 
-S(Q^) = [Z {Q'^) / Z i?(/i^) which reduce in the one-loop approxima- 
tion to Z{Q^) ~ a'^[L] with ly = 7o/(2/3o) being a fractional number. 

All that means that in order to generalize the APT in the "analytization as 
a whole" direction, one needs to construct analytic images of new functions: 
a'^ , L"^, .... This task was performed in the framework of the so-called 
FAPT suggested in [2H]. Now we briefly describe this approach. 

In the one-loop approximation, using recursive relation (|27c|) we can 
obtain explicit expressions for Ai,[L] and 2li^[L]: 

ML\ = j^- ^ (29a) 



sm 



_ — 1) arccos ( L/V^r^ + L'^ 

ML] = ^ , ,,,, . (29b) 

7r(i/- l)(^2 + ^2){-l)/2 



Here F{z, u) is the reduced Lerch transcendental function which is an an- 
alytic function in v. The couplings Au[L] and have very interesting 
properties, which we discussed extensively in our previous papers [25] . 

In Fig. [SI we show in comparison how ^u[L] and Au[L] depend on L 
for fractional values of ly: one more time we observe the same picture 
of the Distorting Mirror when comparing the Minkowski (left panel) and 
Euclidean (right panel) regions. 

To demonstrate the importance of taking into account the FAPT, that 
is using Au[L] and ^u[L] instead of (^i[L])'^ and (2li[L])'', we show in Fig. [9| 
the values of the normalized deviations AuiLju) = l — {%i[L])'^/^^[L] and 
A-^iLjiy) = 1 — (^i[L])'^ /^,y[L] in the Minkowski and Euclidean domains, 
respectively. 



18 



Alexander P. Bakulev, Dmitry V. Shirkov 




2468 10 2468 10 



Figure 9: Left panel: Comparing 2t,y with (Sli)"^ for fractional v — 0.62 (solid 
line), 1.62 (dotted line) and 2.62 (dashed line). Right panel: Comparing Ay with 
(^Aif for the same fractional values of v as in the left panel. 

The construction of the FAPT with a fixed number of quark flavors, 
N is a two-step procedure: we start with the perturbative result [a((5^)] ^ ^ 
generate the spectral density 'Pi,{cr) using Eq. (j26p . and then obtain analytic 

couplings and 2lj^[L] via Eqs. (j27|) . Here Nr is fixed and factorized 

out. We can proceed in the same manner for A'j-dependent quantities: 

[a,{Q^;Nf)Y p,{a;Nf) = p,[L„;Nf] = pM)/b'j ^ A,[L;Nf] and 
2tj,[L; Nf] — here Nf is fixed but not factorized out. 

The global version of the FAPT (2009), which takes into account 
heavy-quark thresholds, is constructed along the same lines but starting 
from global perturbative coupling [a|'°''((5^)] , being a continuous function 

of due to choosing different values of QCD scales Aj, corresponding to 
different values of Nf. We illustrate here the case of only one heavy-quark 
threshold at = m|, corresponding to the transition Nf = 3—7- Nf = 4. 
Then we obtain the discontinuous spectral density 

-^(a) = e{L,< U) pn [L„; 3] + 9 {U < L,) p„ [L, + A4; 4] , (30) 

with L„ = In (a/A2), L/ = In (m)/Kfj and A/ = In (a^/A^) for / = 4, 

which is expressed in terms of the fixed-flavor spectral densities with 3 and 
4 flavors, Pn[La', 3] and Pn[La + A4; 4]; note here that Lg- + A4 = ln((T/A|). 
However, it generates the continuous Minkowskian coupling 

Stf^iL] = e(L<L4) (^%[L; 3] + A43%) + e {L4<L) 2l^[L + A4;4] (31a) 

with A432t^ = 21^ [L4 + A4 ; 4] — 21^ [L4 ; 3] and the analytic Euclidean coupling 

Ar[L]=AAL + X.-A]+ + dL^. (31b) 

J \ + 

—00 
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To demonstrate the magnitude of the threshold corrections, we show in 
Fig. do] the values of the normalized deviations AA„[L] = Af^°^[L] - A^[L+ 
A4; 4] and A2l^ [L] = 21^,'°*^ [L] - %y [L+A4; 4] inthe Euclidean and Minkowski 
domains, respectively (for more details see 



2.4. Electromagnetic pion form factor at NLO 

The scaled hard-scattering amplitude truncated at the next-to-leading or- 
der (NLO) and evaluated at renormalization scale = XrQ^ reads (see 
for details [29]) 



+ 



(ArQ2) 



An 



Crtjf (^'2/;^) +^o4''''^(x,?/;Ar) + 4^^)(x,2/)|(32) 



with shorthand notation {x = 1 — x) 



''H,2 \ Q2 



4°^ ix,y) 



2(3 + ln(xy)) In 



The leading twist-2 pion distribution amplitude (DA) [30] at the normal- 
ization scale is given by [31] 



ip^[X,flp) 



6x{l 



l + Y,a2n{f-ihcli\2x-l) 



n>l 



All nonperturbative information is encapsulated in the Gegenbauer coeffi- 
cients a2n(^|')- To obtain the factorized part of the pion form factor (FF), 
one needs to convolute the pion DA with the hard-scattering amplitude: 



i^J-*(Q^)=<^^(x;^^)C^DriJ 



NLO 



y 
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0.1 , . 1 



10 20 30 40 50 10 20 30 40 50 

[GeV] [GeV'l 

Figure 11: Factorizcd pion FF in the "Naive Analytization" (left panel) and in 
the "Maximal Analytization" (right panel). Solid lines correspond to the scale 
setting = 1 GcV^; dashed lines — to fij^ = Q^, dotted lines — to the BLM 
prescriptions, whereas dash-dotted lines — to the a„-scheme. 



In order to obtain the analytic expression for the pion FF at the NLO the 
so-called "Naive Analytization" was suggested in [32]. It uses the analytic 
image only for coupling itself but not for its powers. In contrast and in 
full accord with the APT ideology the receipt of "Maximal Analytization" 
has recently been proposed in [JU], which uses the analytic image A2 for 
the second power of coupling as well. In Fig. [11] we show the predictions 
for the factorized pion FF in the "Naive" and the "Maximal Analytization" 
approaches. We see that in the "Maximal Analytization" approach the ob- 
tained results are practically insensitive to the renormalization scheme and 
scale-setting choice (already at the NLO level). It is interesting to note here 
that the FAPT approach, used in [23] for analytization of the ln{Q'^ / fi'p)- 
terms in the hard amplitude (j32p . diminishes also the dependence on the 
factorization scale setting in the interval fj.'j^ = 1 — 50 GeV^. 

Using the FAPT it appears to be possible to estimate the NNLO cor- 



Q'-F^{Q'') [GcV^] 1 






2 4 6 


[GcV^] 

8 1 



Figure 12: We show as a narrow dashed-dottcd strip the predictions for the pion 
FF obtained using the improved Gaussian model of the nonlocal QCD vacuum. 
The width of the strip is due to the variation of the Gegenbauer coefBcients 02 and 
a\ in the corresponding shaded bands for the pion DA (indicated by the central 
solid line). Note that this dashed-dotted strip shows the effect of the 0{A-i) 
correction only for the central solid curve of the shaded band. 
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rection to the whole pion FF without very comphc ated calculation of the 
corresponding three- loop triangle spectral density [33]. Here we show in 
Fig. [12] only the results: The NNLO correction is of the order of 3 — 10%. 

2.5. Resummation in the one-loop APT and FAPT 

We consider now the perturbative expansion of a typical physical quantity, 
like the Adler function and the ratio R, in the one-loop APT. Due to limited 
space of our presentation we provide all formulas only for quantities in the 
Minkowski region: 



n[L] 



oo 

E 

n=l 



(33) 



We suggest that there exists a generating function P{t) for the coefficients 
dn = dn/di: 



roo roc 

dn= I P{t) f-^dt with / P{t) dt = l. 
Jo Jo 



(34) 



To shorten our formulae, we use for the integral J^f{t)P{t)dt the following 
notation: ((/(t)))p(t). Then the coefficients dn = di {{t'^~^)) p[t) S'^d, as has 



been shown in [24], we have the exact result for the sum in (|33p 

7^[L] = di((2ti[L-t]))p(i). (35) 

The integral in variable t here has a rigorous meaning ensured by the finite- 
ness of the coupling 2ti[t] < 1 and fast fall-off of the generating function 

In our previous publications |25j , we constructed generalizations of these 
results, first, to the case of the global APT when heavy-quark thresholds 
are taken into account. Then one starts with the series of type (|33p . where 
2t„[L] are substituted by their global analogs 2l^'°'^[L] (note that due to 
different normalizations of global couplings, 21^'°*^ [L] ~ 2l„[L]//3j, the coef- 
ficients dn should also be changed). Then 

t 



7^s'°'=[L] = di^(L<L4)((A42li[t] +21] 



L- 



»P(t) 



+ di0(L>L4)((2l] 



where Ai^^[t] = % 



L4 + A4 



t//34;4 



L+Xa 



21. 



t 



))p{t) ; 



(36) 



t//33;3 



The second generalization has been obtained for the global FAPT. Then 
the starting point is the series of the type Y^^=o dn '^n+uVA ^^^"^ the result 
of summation is a complete analog of Eq. (j36p with the substitutions 



P{t)^Pu{t) 



p 



V X 



dx 



(37) 
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do do^u[L], 2li[L -t] ^ ^i+u[L - t], and A42li[t] A42t_i+^[t]. All 
needed formulas have also been obtained in parallel for the Euclidean case, 
for details see [25], 



2.6. Applications to Higgs boson decay 

Here we analyze the Higgs boson decay to a bb pair. For its width we have 

r(H ^ bb) = Mh Rs{Mfj) (38) 

4v2vr 

with RsiMfj) = ml{Mjj) Rs{Mjj) and Rs{s) is the i?-ratio for the scalar 
correlator, for details see |23^ I34j. In the one-loop FAPT this generates the 
following nonpower expansioiF^: 

[L] = 3 m^, J 21^- [L] + dtY.§i 21^.0 W | ' (39) 

I n>l ) 

where ~ 9-05 it 0.09 GeV^ is the RG-invariant of the one- loop m^(/i^) 



evolution ml{Q'^) = rh?^^^^ a^o(Q^) with vq = 27o/6o(5) = 1.04 and 70 is the 

quark-mass anomalous dimension. This value itt^^i^ was obtained using the 

one-loop relation |35j between the pole 6-quark mass of [36] and the mass 
mfc(mb). 

We take for the generating function P{t) the Lipatov-like model of |25j 
with {c = 2.4, 13 = -0.52} 

It gives a very good prediction for d^ with re = 2, 3, 4, calculated in the 
QCD PT [M]: 7.50, 61.1, and 625 in comparison with 7.42, 62.3, and 620. 
It is worthwhile to remind here the history of calculating the /3-function 
coefficients /3„, in the ip\ scalar field theory. In [37] the resummation proce- 
dure was suggested on the basis of taking into account four-loop results in 
the MOM scheme. The five-loop calculations of the anomalous dimensions 
72 and 74 for this model in the MS scheme was performed in [38] (72) and 
in |39] (74 was calculated numerically with small errors). In this last paper, 

using the Borel-like technique and the four-loop results in the MS scheme 
the five-loop prediction = 1405 ± 80 was made. The calculated re- 

sult /Js = 1420.69 appeared in the range predicted in |37[. The uncertainties 
of numerical calculations were eliminated by Kazakov |40j using the unique- 
ness method for multiloop calculations — he confirmed numerical results. 



Appearance of denominators tt" in association with the coefficients dn is due to d„ 
normalization. 
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After that, in [H], the errors in the previous results for both 72 and 74 
were revealed. As a result, 72 and 74 were changed, but the value of fi^ 
numerically appears to be practically the same! Resume: resummation 
predictions for the MS-scheme are really in very good accord with the 
five-loop results. 

Then we apply the FAPT resummation technique to estimate how 
good is the FAPT in approximating the whole sum TZs[L\ in the range 
L e [11.5,13.7] which corresponds to the range Mh e [60,180] GeV^ 



.Nf=6 

with Aq^o 



189 MeV and Slf'""^ 



0.122. In this range, we have 



(Le = ln(mf/Ai)) 



ns[L] 



L+A5 



71-/35 ■ 



)>ps 



(41) 



with defined via Eqs. (gOD and <^ 

Now we analyze the accuracy of the truncated FAPT expressions 



7^s[L;A^] = 3 m' 



(1) 



V 



n=l 



glob 



[L] 



(42) 



and compare them with the total sum TZ^ [L] in Eq. (j4ip using relative errors 

Aiv[L] = 1 - ■^s[L; N]/ns[L\. In Fig. [T31 we show these errors for iV = 2, 
= 3, and = 4 in the analyzed range of L G [11, 13.8]. We see that 

already TZs[L;2\ gives accuracy of the order of 2.5%, whereas TZs[L;^\ of 
the order of 1%. Looking at Fig. [13] we understand that only in order to 
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Figure 13: Left panel: The relative errors AAr[i], = 2, 3 and 4, of the truncated 
FAPT in comparison with the exact summation result, Eq. (pij) . Right panel: The 
width r^^^j; as a function of the Higgs boson mass Mh in the rcsummed FAPT. 
The width of the shaded strip is due to the overall uncertainties induced by the 
uncertainties of the resummation procedure and the pole mass error-bars. Both 
panels show the results obtained in the one- loop FAPT |25) . 



have the accuracy better than 0.5%, one needs to take into account the 4-th 
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correction. We verified also that the uncertainty due to P(t)-modelhng is 
small < 0.6%, while the on-shell mass uncertainty is of the order of 2%. 
The overall uncertainty then is of the order of 3%, see in the right panel of 
Fig. [131 that is in agreement with the Kataev&Kim estimations |35j . 



3. Conclusions 

• Perturbation approach in Quantum theory (both the nonrelativistic 
and QFT) suffers from divergence of (asymptotic) series in powers of 
small expansion parameter g. The reason is the nonanalyticity (essen- 
tial singularity) at the origin of the g complex plane. 

Due to this, common power expansions have a lower limit of accuracy 
(as shown in Table 1). In particular, this refers to perturbative QCD. 
There, for some low-energy processes (see Table 2) this lower hmit 
exceeds the experimental error. This is the reason that the nonper- 
turbative means are of utmost importance in many physical situations. 

• In Section 1, two of the nonperturbative instruments, the Renormal- 
ization Group and the Dispersion Relation methods, are outlined; de- 
scription of their application to the actual case of Quantum Chromo- 
dynamics then follows. Here, the particular theoretical construction, 
the Analytic Perturbation Theory (APT), was devised on the turn of 
the century. 

• Section 2 begins with the summary of the APT basic elements. We 
also remind that in the APT one has 

— Universal (loop & scheme independent) IR limit; 

— Practical renormalization-scheme independence; 

— Non-power perturbation expansion over a set of particular func- 
tions An{Q'^) , 2l„,(s) instead of the running coupling powers ; 

— Quick loop convergence that improves the situation with the role 
of higher-loop corrections. 

• As a result of quick loop convergence, we show (see Table 3) that 3- 
and 4-loop APT terms contribute to observables less than 5%, i.e., 
below the current level of data errors. Hence, the two-loop, NNLO 
level is practically sufficient. 

• Then we expose the Fractional APT (FAPT) that provides an effective 
tool to apply the APT approach for renormgroup-improved perturba- 
tive amplitudes. By the example of the pion electromagnetic form 
factor we show that the FAPT delivers minimal sensitivity to both 
renormalization and factorization scale setting. 

• In both the APT and FAPT approaches we describe the resumma- 
tion procedures that produce finite resummed answers for perturba- 
tive quantities if one knows the generating functions P{t) for the PT 
coefficients. Using quite simple model generating function Ps{t) for 

Higgs boson decay H ^ bb we conclude that at N^LO we have accu- 
racy of the order of 1% due to the truncation error and of the order 
of 2% due to the RG-invariant mass uncertainty. 
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